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ABSTRACT. We calculate the mean power spectrum of galaxies using 
published power spectra of galaxies and clusters of galaxies. The mean 
power spectrum of cluster samples as well as APM 3-D, IRAS QDOT, and 
SSRS+CfA2 galaxy samples has a relatively sharp maximum at wavcnumbcr 
k = 0.05 ±0.01 h Mpc -1 , followed by an almost exact power-law spectrum of 
index n ?s —1.9 toward smaller scales. The power spectrum found from APM 
2-D galaxy distribution and from LCRS and IRAS 1.2 Jy surveys is flatter 
around the maximum. We show that power spectra of galaxies and matter are 
similar in shape for a wide class of assumptions, the bias parameter of galaxies 
relative to matter is fixed by the fraction of clustered matter associated with 
galaxies; and find b ga [ = 1.3 ± 0.1. We compare the empiric power spectrum 
of matter with analytical power spectra and show that the primordial power 
spectrum has a break in amplitude and a spike if presently available galaxy 
and cluster samples represent the true mass distribution of the Universe. 



1 The power spectrum of 
galaxies and clusters 

There exist a large body of ob- 
servational determinations of the 
power spectrum of galaxies and 
clusters of galaxies. In this talk I 
shall review the determination of 
the mean galaxy power spectrum 
(Einasto etaZ1998a, hereafter E98a). 
The mean galaxy power spectrum 
shall be reduced to that of mat- 
ter (Einasto eiaZ1998b, E98b). This 
semi-empirical matter power spec- 
trum shall be used to determine the 
primordial power spectrum (Einasto 
eta£1998c, E98c). 

Recent determinations of power 
spectra for large galaxy and cluster 
samples are plotted in Figure 1. The 
compilation in based on summary by 
E98a. Spectra are shifted in ampli- 
tude to match the amplitude of the 
power spectrum of APM galaxies on 
wavenumber k = 0.1 h Mpc -1 . The 
APM galaxy spectrum is a recon- 
struction of the 3-dimensional spec- 
trum based on deep 2-dimensional 
distribution of over 2 millions of 
galaxies, thus the cosmic error is 
smaller here than in available 3-D 
surveys. The APM galaxy spectrum 
is also free of redshift distortions. We 
see that after vertical scaling there 
is little scatter between individual 
determinations of power spectra on 
medium and small scales. On large 
scales around the maximum the scat- 



ter is much larger. 

We have formed two mean power 
spectra. One spectrum is based on 
samples having power spectra with 
a high amplitude near the maxi- 
mum. Such samples are Abell-ACO 
and APM cluster surveys, the red- 
shift survey of APM galaxies, and 
the SSRS-CfA2-130 galaxy survey. 
These samples cover large regions in 
space where both high- and medium- 
density regions are well represented, 
thus we use the notation Phd 
for this mean power spectrum (HD 
for high-density). This mean power 
spectrum has a relatively sharp max- 
imum at k = 0.05 ± 0.01 h Mpc -1 , 
followed by an almost exact power- 
law spectrum of index n —1.9 to- 
ward smaller scales. 

The other mean spectrum is based 
on samples which have a power 
spectrum with a shallower turnover; 
such samples are the LCRS survey, 
IRAS QDOT galaxies, and the APM 
2-D sample of galaxies. In LCRS 
and IRAS QDOT surveys medium- 
density regions are well present but 
not regions of highest density (very 
rich superclusters). We use the nota- 
tion Pmd for this mean power spec- 
trum (MD for medium-density). On 
medium and small scales it coincides 
with the previous spectrum, but it 
has a maximum of lower amplitude. 

Presently it is not clear, whether 
the difference between these two 
spectra is real or partly due to some 
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Figure 1. Power spectra of galaxies and clusters of galaxies scaled to match 
the amplitude of the 2-D APM galaxy power spectrum (E98a). Spectra are 
shown as smooth curves and are designated as follows: ACO-E and ACO- 
R - Abell-ACO clusters (Einasto etall997a, Retzlaff etall998); APM-T - 
APM clusters (Tadros etall998); APM-TE - APM galaxies (Tadros & Efs- 
tathiou 1996); APM-P-GB - spectra derived from 2-D distribution of APM 
galaxies (Peacock 1997, Gaztanaga & Baugh 1997); IRAS-P - IRAS galax- 
ies (Peacock 1997); CfA2 - SSRS-CfA2 130 Mpc galaxy survey (da Costa 
etall994); LCRS - LCRS survey (handy etall996); P(k) me an indicates the 
mean power spectrum Phd f or high-density regions; the power spectrum for 
medium- density regions, Pmd, is identified with the spectrum APM-P-GB. 
The mean error of the mean spectrum is 11 %, for individual samples it 
varies between 23 and 48 %■ 



artifacts of data handling. Several 
arguments suggest that there exist 
real differences between power spec- 
tra of various samples. All samples 
which have a high amplitude of the 
power spectrum near the maximum 
are deep fully 3-D samples. In con- 
trast, samples with a shallower power 
spectrum arc not so deep, or do 
not contain very rich superclusters as 
the LCRS sample. In IRAS QDOT 
sample galaxies in rich superclusters 
were removed (Tadros & Efstathiou 
1995). 

On the other hand, some artifacts 
of data reduction or influence of the 
sample selection and/or geometry 
are not excluded. For instance, the 
Las Campanas survey is not a fully 
3-dimcnsional sample, it is made in 
narrow strips which may smooth out 
sharp features of the power spectrum 
near the maximum. A curious fact is 
the reconstruction of the 3-D spec- 
trum from APM 2-D galaxy distri- 
bution which has near the maximum 
a lower amplitude as expected from 



real 3-D observations of galaxies of 
the same sample. This difference is 
not explained yet, it may be due to 
some problems with the reconstruc- 
tion of the 3-D power spectrum from 
2-D data. 

The difference between two mean 
power spectra can be considered as 
the combined result of the cosmic 
scatter and our ignorance of all de- 
tails of the data reduction. 

2 The reduction of galaxy 
power spectrum to matter 

To compare the observed power 
spectrum with theoretical spectra 
the galaxy spectrum must be re- 
duced to matter. 

Differences between the distribu- 
tion of galaxies and matter are due 
to the gravitational character of the 
evolution of the Universe. As shown 
already by Zeldovich (1970), the 
evolution of under- and over-dense 
regions is different. Matter flows 
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Figure 2. The biasing parameter as 
a function of the wavenumber k for 
2-D simulation, determined for all 
galaxies (threshold density po = 1 in 
units of the mean density), galaxies 
in high-density regions (po = 2), and 
galaxies in clusters (po = 5). 
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Figure 3. The fraction of matter 
associated with galaxies, F ga i, as a 
function of time, measured through 
the erg parameter (curved lines). 
Thick straight line shows the rela- 
tion equation (3); open circle notes 
the observed value of erg. 



away from low-density regions to- 
ward high-density ones until it col- 
lapses. In order to form a galaxy or 
a system of galaxies, the mean den- 
sity of matter in a sphere of radius 
r must exceed the mean density by 
a factor of 1.68 (Press & Schechter 
1974), the radius r determines the 
mass of the contracting object. Thus 
in low-density regions (voids) galax- 
ies are absent, but the gravity is 
unable to evacuate voids completely 
— there exists primordial matter in 
voids. Visible matter is concentrated 
together with dark matter in a web 
of galaxy filaments and superclus- 
ters (Zeldovich, Einasto & Shandarin 
1982, Bond, Kofman & Pogosyan 
1996). 

These considerations show that 
model particles can be divided into 
two populations, the unclustered 
population in voids, and the clus- 
tered population in high-density re- 
gions. The last population is associ- 
ated with galaxies including DM ha- 
los of galaxies and clusters of galax- 
ies. To get the clustered population 
one has to exclude the population 
of low-density particles using a cer- 
tain threshold density, po, which di- 
vides the matter into the matter in 
voids and the clustered matter. Hy- 
drodynamical simulations by Cen & 
Ostriker (1992, 1998) show that the 
overall mean density is a good ap- 
proximation to the threshold density. 
In determining the density field we 
use smoothing on scales comparable 



to the size of actual systems of galax- 
ies (~ 1 h^ 1 Mpc). 

Analytical calculations and nu- 
merical simulations show that the ex- 
clusion of matter from low-density 
regions rises the amplitude of the 
power spectrum but not its shape 
(E98b). Power spectra of galaxies 
and matter are related as follows: 

P 9 al(k)=b 2 P m (k), (1) 

where the bias factor b is expressed 
by the fraction of matter in the 
clustered population associated with 
galaxies, 

b = l/F gal . (2) 

Actually the biasing parameter b is a 
slow function of the wavenumber k, 
but in the range of scales of interest 
for the present study it is practically 
constant. The biasing parameter b(k) 
is shown in Figure 2 for a 2-D model, 
it is found by comparing power spec- 
tra of all particles (matter) and clus- 
tered particles associated with galax- 
ies using different threshold density 
po (E98b). 

The fraction of matter in the 
clustered population can be deter- 
mined from numerical simulations of 
the void evacuation. This has been 
done for several models: the stan- 
dard CDM model (SCDM, Hubble 
constant h = 0.5), model with cos- 
mological constant (LCDM, h = 0.7, 
Cl = 0.3), and open model (OCDM, 
h = 0.7, fi = 0.5). In Figure 3 
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Figure 4. The semi- empirical matter power spectra compared with theoreti- 
cal and primordial power spectra for mixed DM models. Left: present power 
spectra; right: primordial spectra. Solid bold line shows the matter power spec- 
trum found for regions including rich superclusters, Ph D(k) ; dashed bold line 
shows the power spectrum of matter Pmd(^) f or medium dense regions in the 
Universe. On small scales observed power spectra are corrected for non-linear 
effects. Model spectra with Qo = 0.9, . . . 0.3 are plotted with solid lines; for 
clarity models with Qo = 1-0 and 0.5 are drawn with dashed lines. Primordial 
power spectra are shown for the power spectrum PijD{k); they are reduced to 
scale-free spectrum, P{k) ~ k. 



we show the increase of the frac- 
tion of matter associated with galax- 
ies for these models (for LCDM 2 
models with different realization are 
given). The epoch of simulation is 
measured in terms of erg, rms den- 
sity fluctuations in a sphere of radius 
r = 8 h Mpc. From model data 
alone it is impossible to determine 
erg which corresponds to the present 
epoch. The value of this parameter is 
known from observations for galax- 
ies. It is found by the integration 
of the observed mean power spec- 
trum of galaxies; we get (crg) 9ai = 
0.89 ± 0.05. This value is equal to erg 
for matter, if all matter were asso- 
ciated with galaxies (F ga i = 1). In 
reality they are different and related 
via an equation similar to (1) since 
cr| is proportional to the amplitude 
of the power spectrum: 

(o"g)m = Fgal(°~s)gal- (3) 

This relationship is plotted in Fig- 
ure 3 by a straight line. The intersec- 
tion of this line with curves F ga i vs. 
erg for models yields values of both 
parameters which correspond to the 
present epoch. We get erg = 0.68 ± 
0.06 for matter; F gal = 0.75 ± 0.05, 
and the biasing parameter of galaxies 
in respect to matter b ga i = 1.3 ± 0.1 
(E98b). 



3 The primordial power spec- 
trum 

The final step in our study is the 
comparison of the power spectrum 
of mass with theoretical power spec- 
tra for models of structure formation. 
We use models with cold dark mat- 
ter (CDM) and a mixture of cold and 
hot dark matter (MDM) in spatially 
fiat models. We derive the transfer 
functions for a set of cosmological pa- 
rameters. 

Models arc normalized on large 
scale by four-year COBE normaliza- 
tion (Bunn & White 1997); the den- 
sity of matter in baryons is taken 
Qf, = 0.04 (in units of the critical 
density); and the Hubble parame- 
ter h = 0.6. The cosmological con- 
stant was varied between Q/± = 
and Qa = 0.8. In mixed dark mat- 
ter models the density of hot DM was 
fixed, Q v = 0.1, cold DM density was 
chosen to get a spatially flat model. 

Analytical power spectra for 
MDM models arc plotted in the left 
panel of Figure 4 together with the 
semi-empirical matter power spectra, 
P H D(k) and P MD (k). MDM models 
fit the semi-empirical matter power 
spectrum better than other models. 
In the right panel we show the initial 
power spectrum, 

Pi„it(k) = P(k)/T 2 (k), (4) 
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compared with the scale-free primor- 
dial power spectrum, Po(fc) ~ k; here 
T(k) is the transfer function. The ini- 
tial power spectrum is plotted only 
for spectrum PuD{k) which corre- 
sponds to high-density regions. 

The main feature of primordial 
power spectra is the presence of a 
spike at the same wavenumber as 
that of the maximum of the ob- 
served power spectrum. Primordial 
spectra also have a break, i.e. their 
amplitudes on small scales are dif- 
ferent from the amplitudes on large 
scales. The shape of the primordial 
spectrum varies with the cosmologi- 
cal constant. The primordial power 
spectrum which is found from the 
shallower spectrum Pmd(') has no 
sharp peak but the break is similar 
to that for the spectrum Prd- 

The main conclusion from the 
present analysis is that it is impos- 
sible to avoid a break and/or spike 
in the primordial power spectrum, 
if presently available cluster and 
galaxy data represent a fair sample 
of the Universe. Clusters of galaxies 
cover a much larger volume in space 
than galaxies, thus cluster samples 
are presently best candidates for the 
fair sample. However, this conclu- 
sion is tentative. New very deep sur- 
veys of galaxies now in progress will 
specify the power spectrum on large 
scales more exactly and yield a bet- 
ter estimate of the primordial power 
spectrum. Presently we can say that 
the possibility of a broken scale and 
peaked initial power spectrum has to 
be taken seriously. 
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